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^ ^ Abstract. For a sample of nine Galactic globular clusters we measured the inner kinematic 

^ profiles with integral-field spectroscopy that we combined with existing outer kinematic 

measurements and HST luminosity profiles. With this information we are able to detect 
the crucial rise in the velocity-dispersion profile which indicates the presence of a central 
^ \ black hole. In addition, N-body simulations compared to our data will give us a deeper 

, insight in the properties of clusters with black holes and stronger selection criteria for further 

^>«^ ' studies. For the first time, we obtain a homogeneous sample of globular cluster integral- 

' field spectroscopy which allows a direct comparison between clusters with and without an 

, intermediate-mass black hole. 
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1. Introduction X-ray and radio emissions, the central kine- 
matics in globular clusters can reveal possi- 

d' T 1- 11111 ^T»iT-.rT ble IMBHs. However, resolving the gravita- 

M ' Intermediate-mass black holes (IMBHs, M, ~ i . £■ a r i i-i Tivjr.TT 

Cn i , ^9 X 1 1 1 • tional sphere or influence tor plausible IMBH 

- - 10 - 10 Mo) have drawn the attention or as- c i i 11.1 

J. , , , , , , masses (1 - 2 tor large Galactic globular 
tronomers tor more than a decade and would 1 , n • 1 J• 
, , ,. , , ^, ., , clusters) requires velocity-dispersion measure- 
shed light onto the mystery or the rapid growth . 1 u- u 1 1 ■ 
^ - iiiii , . ^ , ments at a high spatial resolution. Today, with 



of supermassive black holes by acting as seeds 
in the early stage of galaxy formation. Besides 



the existence of the Hubble Space Telescope 
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(HST) and high spatial resolution ground 
based integral-field spectrographs, the search 
for IMBHs is possible. 

A strong motivation for searching for 
intermediate-mass black holes is the observed 
relation between the black-hole mass and the 
velocity dispersion of its host galaxy (e.g. 
Ferrarese & Merritt 2000; Gebhardt et al. 
2000a; Giiltekin et al. 2009). Exploring this re- 
lation in the lower mass range, where we find 
intermediate-mass black holes, will give im- 
portant information about the origin and vali- 
dation of this correlation. 

The formation of intermediate-mass black 
holes can occur by the direct collapse of very 
massive first generation stars (Population HI 
stars, Madau & Rees 2001), or runaway merg- 
ing in dense young star clusters (e.g. Portegies 
Zwart et al. 2004). This makes globular clus- 
ters excellent environments for intermediate- 
mass black holes. The currently best candi- 
dates for hosting an intermediate-mass black 
hole at their centers are the globular clusters 
(x> Centauri (NGC 5139, Noyola et al. 2008, 
2010; van der Marel & Anderson 2010; Jalali 
et al. 2012), Gl in M31 (Gebhardt et al. 2005), 
and NGC 6388 (Liitzgendorf et al. 2011). 
AH of these very massive globular clusters 
show kinematic signatures of a central dark 
mass in their velocity-dispersion profiles. For 
(jj Centauri, van der Marel & Anderson (2010) 
investigated a large dataset of HST proper 
motions and found less evidence for a cen- 
tral IMBH than proposed by Noyola et al. 
(2008). Using a new kinematic center, how- 
ever, Noyola et al. (2010) and Jalali et al. 
(2012) confirmed the signature of a dark mass 
in the center of oj Centauri and proposed a 
-40 000 Mo IMBH. 

Further evidence for the existence of 
IMBHs is the discovery of ultra luminous X- 
ray sources at non-nuclear locations in star- 
burst galaxies (e.g. Matsumoto et al. 2001; 
Fabbiano et al. 2001). The brightest of these 
compact objects (with L ~ 10"^' ergs"') imply 
masses larger than IO^Mq assuming accretion 
at the Eddington Umit. 



2. Observations 

Our sample consists of nine Galactic globu- 
lar clusters, each of them a good candidate for 
hosting an IMBHs at its center due its high 
mass and central cusp in its fight profile. 

These clusters were observed with the 
GIRAFFE spectrograph of the FLAMES 
(Fiber Large Array Multi Element 
Spectrograph) instrument at the Very Large 
Telescope (VLT) using the ARGUS mode 
(Large Integral Field Unit). The velocity- 
dispersion profile was obtained by combining 
the spectra in radial bins centered around the 
adopted photometric center and measuring the 
broadening of the lines using a non paramet- 
ric line-of-sight-velocity-distribution fitting 
algorithm (for more detailed description see 
Lutzgendorf et al. 2011, 2012a; Lutzgendorf 
et al. 2012c, and Feldmeier et al. 2013, in 
prep.). For larger radii, the kinematic profiles 
were completed with radial velocity from 
the literature, if existent. In addition to the 
spectroscopic data, HST photometry was used 
to obtain the star catalogs, the photometric 
center of the cluster and its surface brightness 
profile. For each cluster both photometry 
and spectroscopy were combined in order 
to apply analytic Jeans models to the data. 
The surface brightness profile was used to 
obtain a model velocity-dispersion profile 
which was fit to the data by applying diff'erent 
black-hole masses and M/Ly profiles. The 
final black-hole masses were obtained from a 
)^ fit to the kinematic data. 



3. N-body simulations 

Furthermore, we run N-body simulations 
based on the GPU (Graphic Processing Unit) - 
enabled version of the coUisional N-body code 
NBODY6 (Aarseth 1999). The code treats bi- 
nary interactions, stellar evolution, and exter- 
nal tidal fields in a highly sophisticated man- 
ner which allows detailed simulations of glob- 
ular clusters with and without central IMBHs 
and their observable effects. The goals of these 
simulations can be summarized in three major 
points: 
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Fig. 1. Velocity maps of the nine Galactic globular clusters of our sample. 



1 ) Reproduce Observations: Our observed 
clusters will be compared to a grid of scaled 
N-body simulations in order to find the best fit- 
ting initial parameters and constrain the black- 
hole mass, M/Ly profiles and mass functions. 
In contrast to Jeans models, N-body simula- 
tions can consider the effect of bright stars on 
the velocity measurements as well as strong 
mass segregation of stellar remnants in the cen- 
ter 2) Test our Analysis: With the output of 
the N-body simulations we are creating mock 
datasets of HST images and IFU datacubes in 



order to test our analyzing methods in photom- 
etry and spectroscopy. This will give new in- 
sights on how reliable the method of integrated 
light is in terms of recovering the inner dynam- 
ics of dense stellar systems. 3) Study cluster 
properties: We are analyzing cluster simula- 
tions with different black-hole retention frac- 
tions, IMBH masses and binary fractions in 
tidal fields. With these models we will be able 
to put constraints on our preliminary black- 
hole masses, verify the measurements and pre- 
dict new observables for finding IMBHs. 
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Fig. 2. M. - a relations of IMBHs and SMBHs 
in comparison. The slope of the bestfit to the GCs 
(green line) is by a factor of two smaller than the 
slope of the SMBHs in galaxies (blue line). 



4. Results and Conclusions 

We have investigated the presence of IMBHs 
in nine Galactic globular clusters using a 
combination of HST photometry and integral- 
field spectroscopy. Comparing the velocity- 
dispersion profiles in the central region of 
each GCs with analytical Jeans models yielded 
four candidates which show an IMBH signa- 
ture in their kinematics (see Liitzgendorf et al. 
2011, 2012a; Lutzgendorf et al. 2012c, and 
Feldmeier et al. 2013, in prep.). Using our re- 
sults and IMBH measurements from the litera- 
ture, we study the M, - cr relation for IMBHs 
and compare it to supermassive black holes in 
galaxies. Figure 2 (Liitzgendorf et al. 2013a, 
in prep.) shows this comparison and indicates 
a lower slope for the coiTelation at the lower 
mass end compared to the slope obtained for 
the SMBHs by McConnell et al. (2011). The 
reason for this behavior could be explained by 
high mass loss of the globular clusters due to 
tidal stripping. This could lower the velocity 
dispersion of the system and moves the globu- 
lar clusters away from the M, - cr relation. 

In addition to the observational work we 
performed N-body simulations in order to ad- 
dress open question and degeneracy. The simu- 



lations include stellar evolution, binary forma- 
tion and an external tidal field while the IMBH 
mass, the black-hole retention fraction and the 
primordial binary fraction are varied. With the 
outcome of these simulations we will be able 
to reproduce our observations, test our analy- 
sis methods, and get a deeper understanding of 
the observable effects of IMBHs to their envi- 
ronment (Liitzgendorf et al. 2013b, in prep.). 
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